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Chapter  1:  INTRODUCTION 


Millimeter-wave  systems  have  received  increased  interest  in  recent  years,  and  serve  a  broad  spec¬ 
trum  of  applications.  Many  of  these  systems  employ  receiver  and  transmitter  sub-systems.  The  use  of 
integrated  circuits  at  millimeter- wave  frequencies  is  essential  to  meet  performance,  cost,  size,  and  other 
constraints.  The  integration  of  millimeter-wave  transmitters  and  receivers  is  of  primary  importance. 
It  is  desirable  to  place  as  many  components  as  possible;  for  example,  amplifiers,  filters,  switches,  anten¬ 
nas;  on  one  substrate  or  waveguide  medium.  Millimeter-wave  integrated  circuits  are  receiving  consid¬ 
erable  exposure  [l],  [2]. 

This  thesis  considers  several  important  elements  for  use  in  millimeter-wave  integrated  circuits. 
Two  different  wave-guiding  media  are  considered:  fin-line,  which  is  a  member  of  the  metallized  sub¬ 
strate  family,  and  image  guide,  which  is  a  member  of  the  rectangular  dielectric  waveguide  family. 

Filter-type  structures,  which  consist  of  a  series  of  discontinuities,  are  studied.  These  structures 
are  used  in  matching  elements  for  active  devices  and  antennas,  radiating  elements,  and  as  filters.  Some 
of  these  components  have  been  studied  [3H7l 

The  solution  of  the  discontinuity  problem  is  therefore  of  particular  importance.  A  suitable  char¬ 
acterization  for  a  single  discontinuity  permits  the  analysis  of  some  series  of  discontinuities  via  general¬ 
ized  scattering  and  transmission  matrices. 

The  uniform  fin-line  structure  is  analyzed  in  Chapter  2.  This  analysis  gives  propagation  con¬ 
stants,  characteristic  impedances  and  mode  functions. 

Chapters  3  and  4  deal  with  two  different  approaches  for  solving  the  fin-line  discontinuity  prob¬ 
lem.  A  moment  method  is  used  to  study  a  single  and  an  infinite  periodic  array  of  discontinuities  in 
Chapter  3.  An  iterative  method  and  a  variational  approach  are  used  to  solve  the  single  fin-line  discon¬ 
tinuity  problem  in  Chapter  4. 

An  image  guide  grating  filter  is  studied  in  Chapter  5.  Theoretical  and  experimental  curves  for  the 
insertion  loss  charac  censtic  are  given. 


CHAPTER  2:  UNIFORM  FIN-LINE  ANALYSIS 


2.1.  Introduction 

The  unilateral  fin-line  structure  shown  in  Fig.  2.1.  is  studied  in  this  thesis.  It  consists  of  a  metal 
cladding  on  one  side  of  a  dielectric  substrate,  which  is  mounted  in  a  shielded  enclosure.  The  metalli¬ 
zation  will  be  assumed  infinitesimally  thin.  A  generalized  guide  may  consist  of  some  combination  of 
dielectrics  and  metallized  surfaces. 

A  number  of  recent  papers  have  used  the  spectral  domain  formulation,  in  conjunction  with  a 
moment  method,  to  analyze  the  fin-line  [8H12].  In  particular,  Galerkin’s  method  has  been  shown  to 
be  satisfactory. 

Several  of  these  authors  have  analyzed  the  fin-line  structure  by  expressing  the  fields  in  terms  of 
scalar  potential  functions  in  each  transverse  region,  and  then  applying  the  appropriate  boundary  con¬ 
ditions  to  obtain  the  required  coupled  equations  [8l  [91  [ill  This  approach  is  satisfactory  when  the 
structure  can  be  analyzed  as  a  two-region  problem,  for  example,  symmetric  bilateral  fin-line.  How¬ 
ever,  the  unilateral  fin-line  is  a  three-region  problem,  which  results  in  prohibitive  algebra.  For  a 
larger  number  of  regions  the  problem  is  compounded. 

An  alternative  approach  is  to  use  the  transverse  resonance  technique  in  the  spectral  domain  [lOl 
[121  This  method  greatly  simplifies  the  formulation,  and  will  be  utilized  in  the  work  presented  here. 

The  uniform  fin-line  problem  is  formulated  in  Section  22.  Details  of  the  formulation  procedure 
and  the  numerical  methods  used  in  the  solution,  namely,  the  moment  method,  are  given.  It  is  neces¬ 
sary  to  obtain  data  for  the  propagation  constants,  characteristic  impedances  and  the  mode  functions. 
Theoretical  and  experimental  results  are  given  in  Section  2.3. 

22.  Spectral  Domain  Formulation  for  Fin-Line 

Waveguides  such  as  fin-line  and  microstrip  are  most  simply  formulated  in  the  spectral  domain. 
This  involves  using  the  Fourier  series  representation  of  all  field  quantities  in  the  transverse  direction. 


a 


4 


parallel  to  the  metallization.  This  is  the  x-direction  for  the  geometry  of  Fig.  2.1.  The  transform  pair 
for  quantity  A  is  defined  as 
2* 

A(kxty)=  —  /  A(x,y)eJt“*  dx  (2.1a) 

2*  o 

A  (x,y  )*'*'*  =  S  Att^y)*”''*’* (2.1b) 

/I  *— OO  t 

where  *,  n  =  0,  ±  1,  ±  2- 

The  continuous  Fourier  Transform  may  be  used  to  find  these  quantities  at  the  discrete  kx  values. 

Without  the  metallization  on  the  dielectric  surface,  TE-to-y  and  TM-to-y  modes  propagate.  The 
fin-line  modes  can  be  considered  as  superpositions  of  these  modes. 

An  analogy  is  now  drawn  with  a  plane  wave  incident  on  a  dielectric  boundary  at  some  angle. 
If  the  fields  are  decomposed  into  components  parallel  and  perpendicular  to  the  propagation  direction, 
no  coupling  between  these  components  will  occur  at  the  interface.  The  resulting  equations  have  been 
decoupled. 

A  similar  concept  may  be  employed  in  the  fin-line  case.  Equation  (2.1b)  represents  a  superposi¬ 
tion  of  plane  waves  propagating  in  the  (kx x  +  kzz)  direction,  which  will  be  called  the  u-direction. 
Consider  the  orthogonal  (u,v)  coordinate  system  in  the  (x^)  plane,  as  shown  in  Fig.  2.2(a).  The  fol¬ 
lowing  coordinate  transformation  relationships  hold  for  quantity  a 

au  =  az  Nz  +  ax  Nx  (22a) 

oLy  =  —a.  Nx  +  a,  Nz 

az-auNx+avNz  (22b) 
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Figure  2.2.  (a)  (x^)  and  (u,v)  coordinate  system,  (b)  Equivalent  transverse  resonance  circuits  for 

analyzing  the  fin-line  structure. 


The  field  quantities  expressed  in  the  (u,v)  coordinate  system  allow  a  decoupled  transverse  reso¬ 
nance  analysis.  The  two  equivalent  transmission  lines  for  TM  to  y  and  u  (E fJEuJ2  „)  and  TE  to  y 
and  u  (H ,  JEV  JX  u )  are  shown  in  Fig.  2.2(b).  The  wave  impedances  in  the  y-direction  are 


7  _ 

Eu.  _  >, 

7  — 

—E  „  ;  oj  n 

t-TEi  ~ 

y, 

where 

y,  =  V*,2  +  kz2 

(2.4) 


Spectral  domain  relationships  between  the  fin  current  densities  and  the  slot  fields  may  be 
formed  simply  from  the  transmission  line  equivalent  circuits.  The  relationships  are 

Ju  (*,2,  +  2d)  =  Ym  Eu  ikxXx  +  2d) 

/,  (kzXx  +  2d)  =  Yt  Ev  (kx  Xx  +  2d) 

1  .  Zm  i  +  Ztm2  tanh(y^  i)tanh(2yi<f ) 


(23a) 

(23b) 


where  Ym  = 


Yt 


Z-tmz  tanh(-y2l2)  Ztmi  ^tm 2  tanh(y2l  j)  +  Ztmi  tanh(2'>'  \d )) 
1  ZTE 1  +  ZTE  2  tanh(y2z  1)tanh(2y1<l ) 


Zjmi  tanh(y2l  2)  ZEE  t  (ZEE  2  tanh(y2l  j)  +  ZEE  j  tanh(2y  yd )) 

Transformation  of  Eq.  (23)  back  into  (x^)  coordinates  results  in  two  coupled  equations,  which 


are 


Jz  (kxXz  +  2d)  =  (N*Ym  +  N/Yt)  Ex(kxXy  +  2d)  +  NXNZ  (,Ym  -  Y,  )Ez(kx  X ,  +  2d)  (2.6a) 
/,(*,;,  +  2d  )  =  Nx  Nz  (Ym  -  Yt  )EX  (kx  X  i+2d  )  +  (N*Ym  +  N*Yt  )EZ  ikxiy  + 2d).  (2.6b) 

Equation  (2.6)  is  a  spectral  domain  relationship  between  the  fin  currents  and  the  slot  fields.  This 


may  be  written  in  the  form 


where  the  Y-matm  is  a  type  of  dyadic  Green’s  function.  Only  minor  alterations  to  the  transmission 
line  equivalent  circuits  are  required  to  analyze  an  arbitrary  number  of  dielectric  layers  and  metal¬ 
lized  surfaces. 

The  moment  method  may  be  used  to  obtain  a  satisfactory  numerical  solution  to  Eq.  (2.6).  A 
matrix  equation  may  be  obtained  by  the  application  of  Galerkin’s  method  (testing  functions  being  the 
same  as  the  basis  functions)  to  Eq.  (2.6).  The  slot  fields  are  expanded  in  terms  of  a  suitable  set  of  basis 
functions,  which  satisfy  the  edge  condition,  and  are  readily  Fourier  transformable.  The  representa¬ 
tion  with  P  basis  functions  is 


£,(*)=  £  ap  $,(x) 

pm  i 

£,(*)=  £  \  TJ^x). 


(2.8a) 


(2.8b) 


A  suitable  set  of  functions  is  an  orthogonal  polynomial  set,  modified  by  the  edge  condition.  The 
functions  used  are 


£,(*)  = -  w  — ■  p  =  U3_ 

L  2(x-s)  1 2 


w  ^  .  ,  w 

S - —  <X  <5  +-— 

2  2 


sin  [£—(x  —s  — “)] 

t),(x)  = -  W -  -.<?  =  U4- 

/ 2(x  —s )  2 


£f(x)  =  0 
T),(x  )  =  0 


i  otherwise. 
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Frequently,  £0(.r )  =  Px  (s  ,w )  is  used.  These  functions  are  shown  in  Fig.  2.3  and  Fig.  2.4. 

For  a  centrally  located  slot,  only  even  £  and  odd  7)  terms  (with  respect  to  the  center  of  the  slot) 
are  required,  i.e,  p-0,2,4,_  <7«=2,4,6_  . 

The  matrix  equation  is  formed  using  the  inner  product. 

<*,(*,)//*,)>■  £  x;  (kx  )Y/kz )  (2.io) 

n  =*— oo 

where  kx  = 

2 b 

Of  course,  the  summation  must  be  truncated  to  a  finite  number  of  terms.  The  resulting  homogeneous 
equation  is  of  the  form 

[A]  [x  j  =  0.  (2.11) 


For  P  Ex  basis  functions  and  Q  E.  basis  functions,  A  is  a  (P-t-Q,  P+Q)  matrix  and  x  a  (P+Q,l)  vector 
containing  the  unknown  basis  function  coefficients.  The  eigenvalue  solutions  for  k,  are  found  by 
solving  det  (A)  -  0.  The  coefficient  matrix  x  may  then  be  found.  P-t-Q- 1  unknowns  may  be  found  in 
terms  of  one  unknown.  Thus  the  approximate  Ex  and  E,  eigenfunctions  in  the  slot  can  be  found. 
The  fields  over  the  complete  waveguide  cross-section,  and  hence  the  mode  functions,  may  be  found  by 
returning  to  the  equivalent  transverse  resonance  circuits. 

There  is  no  unique  definition  for  the  characteristic  impedance  for  non-TEM  modes.  A  suitable 
definition  as  discussed  in  [10],  [13]  is 


(2.12) 


The  mode  power  is  given  by 


baa  I fi  function  amplitude 


P  =  Re  \f  f  Ex(x,y)H‘(s,y)  —  Ey(x,y  )H *(x,y )  dxdy  . 


(2.14) 


Using  Parse vaJ’s  theorem,  this  may  be  expressed  as 


2a 

P  =  -j-  Re  E  /  U?* (*x  , y ,y)~Ey(kx  ,y )H‘X (kx  , y )]dy 

n  a  flft  A 


(2.15) 


ExMyJZ,  and  Hx  may  be  found  in  the  (u,v)  coordinate  system  as  voltages  and  currents  to  the 
equivalent  transmission  lines,  giving 


Ex(kx  ,y )H'y(kx  ,y)=  £-1  Ev(kx  ,y)l2+  Eu(kx  .y )E‘v(kx  ,y ) 
£,(*,  ,y)if;a,  ,y)  =  tffB’(*x  ,y)N,  +  if ;(*,  ,y)Nz). 


(2.16a) 


(2.16b) 


The  y-integrals  may  be  performed  analytically.  This  approach  simplifies  the  computation  of  the 
mode  power. 

The  fields  for  the  nth  mode  may  be  expressed  as 


E(x ,y  ,z )  =  e„  (x  ,y  )e  *a*  +  eat(x,y )e 


H(x,y,z)  =  hn(x,y)e~Jkt*z  +  hB,(x,y'k 


2 


(2.17a) 

(2.17b) 


where  e„  and  h„  are  transverse  vector  mode  functions,  and  em  and  hm  are  longitudinal  vector  mode 
functions.  These  mode  functions,  when  normalized,  satisfy  the  following  orthogonality  relation  [14] 

f  f  «*  xhm  zds  =  Snm 

S 

where  Snm  is  the  Kronecker  data  which  is  defined  by 


_  1  n  =  m 

0  n  m 


(.2.1 8) 
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23.  Results  and  Discussion 

For  the  purposes  of  this  study,  a  centrally  located  slot  will  be  considered.  The  techniques 
presented  can  be  used  to  analyze  an  arbitrary  number  of  slots  located  anywhere  in  the  shielded  enclo¬ 
sure. 

Satisfactory  approximate  results  for  the  propagation  constant  can  be  obtained  by  using 
Et  -  Px  (s  ,w )  and  E2  -  0.  Upon  investigation  of  the  convergence  of  the  propagation  constant  solu¬ 
tions,  two  basis  functions  from  Eq.  (2.8)  (  p“l,3,  q-2,4)  give  sufficient  accuracy.  A  "rule  of  thumb" 
for  the  number  of  spectral  terms  required  in  the  inner-product  of  Eq.  (2.9)  is 

"H  X"ikx)Y j(kx)  ,  N2— .  (2.19) 

n  »-.V  /2  W 

This  indicates  that  more  spectral  terms  are  required  for  a  larger  waveguide  width  to  slot  width  ratio. 
This  number  of  basis  functions  gives  good  results  without  excessive  computation  time.  Eigenvalue 
solutions  to  Eq.  (2.10)  can  be  found  by  a  course  iteration  to  find  approximate  solutions,  followed  by 
iterations  using  Newton’s  method. 

The  effective  dielectric  -  kf/kf)  constant  results  agree  with  Schmidt  et  al.  [10].  A  com¬ 
parison  of  the  computed  normalized  wavelength  with  Knorr  and  Shayda  [8]  is  given  in  Fig.  23. 
Knorr  and  Shayda  use  a  pulse  approximation  for  Ex  in  the  slot. 

Experimental  propagation  constants  at  E-band  (76-82  GHz)  have  been  obtained  [15].  Figures  2.6 
and  2.7  give  comparisons  between  the  theory  presented  here  and  this  experimental  data.  The  agree¬ 
ment  is  quite  good.  Additional  propagation  data,  for  propagating  and  evanescent  modes,  which  will 
be  used  in  subsequent  chapters,  is  presented  in  Figs.  2.8,  2.9  and  2.10. 

The  characteristic  impedance,  as  defined  by  Eq.  (2.12),  is  plotted  for  several  cases  in  Figs.  2.11 
and  2.12.  A  pulse  representation  of  Ex  is  used  here.  Note  that  the  impedance  increases  with  increas¬ 
ing  slot  width. 

The  nature  of  the  Green’s  function  matrix  Y  is  of  importance  when  computing  infinite  summa¬ 
tions  and  integrals  in  the  spectral  domain.  The  asymptotic  form  of  Y  is  particularly  important.  For 
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Figure  2.6.  Comparison  of  theoretical  and  experimental  results  for  the  effective  dielectric  constant 
of  fin-line  in  a  WR-12  shield. 
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2.7.  Comparison  of  theoretical  and  experimental  results  for  the  effective  dielectric  constant 
of  fin-line  in  a  WR-28  shield. 
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Figure  2.9.  Dispersion  curves  for  three  even  modes  ( Ex  even  about  x  =$ )  in  unilateral  fin-line 
with  a  WR-28  shield.  / ,  -  l  2  *  3.429  mm,  2d  -  0.254  mm,  €r  ”  2.22,  2b  -  2.556  mm,  w 
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and  for  large  kx 


d-ysgn(kx) 


d  ysgn  (kx ) 

<*4!*,  l+ds 


where  the  c's  and  ds  are  constant  coefficients. 


The  transverse  vector  mode  functions  are 

e(x  ,y )  =  xex  (x  ,y )  +  ye,(x  ,y  ) 
h(x,y)  =  xMx.y)  +  yhy(x,y) 


(221) 


(222a) 

(2.22b) 


These  field  components  are  plotted  for  a  typical  case  in  Figures  2.13  through  2.16.  Notice  that  the 
fields  are  confined  to  the  region  around  the  slot.  The  functions  represent  a  single  propagating  mode. 

The  solutions  for  the  propagation  constant  are  very  stable  when  a  relatively  small  number  of 
basis  functions  and  spectral  terms  are  used,  (i^,  a  small  number  of  spectral  terms  are  required  for 
convergence  in  the  propagation  constant.)  However,  at  least  an  order  of  magnitude  more  spectral 
terms  are  required  for  convergence  of  the  basis  function  coefficients.  Also,  an  increase  in  the  number 
of  basis  functions  necessitates  an  increase  in  the  number  of  spectral  terms. 

It  is  therefore  quite  difficult  to  find  a  satisfactory  representation  for  the  propagating  and  evanes¬ 
cent  mode  functions.  The  analysis  of  a  number  of  fin-line  problems  requires  the  generation  of  a  suit¬ 
able  set  of  orthogonal  modes. 
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CHAPTER  3:  MOMENT  METHOD  ANALYSIS  OF  FIN-LINE  DISCONTINUITIES 
3.1.  The  Fin-Line  Discontinuity  Problem 

Unilateral  fin-line  will  be  considered  in  the  discontinuity  analysis.  The  word  "discontinuity" 
will  mean  a  perturbation  in  the  slot  width.  A  step  discontinuity  is  shown  in  Fig.  3.1a.  The  solution 
of  such  a  problem  is  essential  for  design  and  analysis  of  many  millimeter-wave  components. 

Mode  matching  has  been  used,  with  some  success,  to  analyze  fin-line  discontinuities  [4],  [7],  [l6l 
[17].  An  analysis  of  an  infinite  periodic  array  of  stubs  has  been  published  recently  [18].  This 
approach  utilizes  a  moment  method. 

This  chapter  considers  moment  method  solutions  which  are  applicable  to  a  variety  of  fin-line 
-discontinuities,  including  the  single  discontinuity  of  Fig.  3.1a,  and  the  periodic  array  of  Fig.  3.1b. 

The  formulation  employed  in  this  chapter  reduces  a  three-dimensional  scattering  problem  to  a 
two-dimensional  spectral  relationship  in  the  plane  of  the  fin.  This  technique  employs  the  known 
spectral  Green’s  function,  T .  A  spectral  domain  formulation,  which  is  an  extension  of  the  two- 
dimensional  analysis  of  Chapter  2,  is  employed.  Section  3.2  deals  with  the  single  discontinuity  prob¬ 
lem  and  Section  3  J  considers  the  infinite  periodic  problem.  Conclusions  on  these  techniques  are 
drawn  in  Section  3.4. 

32.  Theoretical  Analysis  for  the  Single  Discontinuity 

The  discontinuity  of  Fig.  3.2  will  be  considered  for  this  analysis.  The  technique  is  quite  general, 
and  may  be  applied  to  other  geometries.  A  series  of  discontinuities,  such  as  Fig.  3.2,  is  of  interest  in 


filtering  and  matching  applications. 

A  spectral  relationship  in  the  "fin-plane"  may  be  written  as  a  general  extension  of  Eq.  (2.6) 


0 


<* 


for  single  in-line  discontinuity. 


Y(kzJtz)E(kzJcz)  =  J(kzJtz) 
where 


(3.1) 


Y(kzJcz)  = 


Ya 


The  spectral  variable  kz  is  again  discrete;  kz  =  2nir/2b,  and  kz  is  a  continuous  spectral  variable  for 
the  single  discontinuity  problem.  Y  is  as  for  Eq.  (2.7),  with  kz  a  spectral  variable  rather  than  an 
eigenvalue. 

A  matrix  equation  can  be  obtained  by  applying  the  moment  method  to  Eq.  (3.1)  with  a  two- 
dimensional  inner  product  defined  by 


<X,(kx*z\Y j{Xz*z)>  *  £  f  X;(kzJct)Yj(kzJcz)dkz.  (32) 

A  single  propagating  mode  will  be  considered  in  this  analysis.  The  desired  information  is  the 
scattering  parameters  of  the  junction.  Consider  a  three  region  problem,  as  shown  in  Fig.  3.2.  The 
fields  in  the  A-  and  C-  regions  are  assumed  known,  except  for  complex  constants.  In  Fig.  32,  a  is 
sufficiently  large  so  that  the  evanescent  modes  are  negligible  in  regions  A  and  C  The  field  distribu¬ 
tions  in  regions  A  and  C  are  found  from  a  uniform  slot  analysis  and  are  themselves  a  superposition  of 
basis  functions.  The  slot  fields  in  the  B-region  can  be  expressed  as  a  superposition  of  suitable  basis 


functions.  We  can  write 


30 


e~ea  +  eb  +  ec 


where 


Q  =  A£C 


Eq  —  Eq  +  Eq ". 


(3.3) 


(3.4) 


The  (+)  signifies  forward  traveling  waves,  and  the  (-)  signifies  negative  traveling  waves.  Expressing 
the  B-region  fields  in  terms  of  a  sum  of  basis  functions, 

Exi  *  £  ^  &  (x,z) 

1*1 

Erf  =  £  bjT)j  (xj).  (3 JS) 

The  fin  current,  ~J ,  and  the  slot  field,  E  ,  are  orthogonal,  and  by  Parseval’s  Theorem,  J  and  E  are 
orthogonal.  Therefore,  the  right-hand  side  of  Eq.  (3.1)  will  not  contribute  in  the  matrix  equation.  Of 
course,  the  matrix  equation  will  have  a  right-band  side  which  is  due  to  the  known  incident  field  on 
the  junction  (£/).  Applying  a  moment  method  to  Eq.  (3.1), 


<in?~Ez  >  +  <(mX^Ez  >  =  0  (3.6) 

<*, Ex  >  +  <+nXuEt>=  0. 

This  equation  is  of  the  form  Ax  -  b,  where  x  is  the  unknown  containing  the  coefficients  of  the 
basis  functions,  and  b  is  due  to  the  incident  field  which  is  known. 

The  reflection  and  transmission  coefficients  at  z  ~  ±  a  can  be  related  to  the  unknown  basis 
function  coefficients  in  the  B-region  by  integral  relationships  on  these  boundaries.  Therefore,  there  are 
I  +•  J  unknowns  to  be  solved  for. 

The  z-dependencies  for  the  incident  (/,),  reflected  ( / r ),  and  transmitted  (  / 1 )  fields  are 


/,(*)»  [l -«_,(*  )]e-'** 

/,(*)«[  l-u_,(r)]e"^,z  (3.7) 

/,tz)  *  . 


The  Fourier  Transforms  of  these,  upon  application  of  the  Weiner-Hopf  technique  (where  a  finite  loss 
is  introduced)  [19],  become 


frOc,) 


f '.(*,) 


-1 

jOtz-^y 


(3.8) 


It  should  be  noted  that  to  reproduce  the  uniform  solution  from  the  formulation  here,  it  is  necessary  to 
consider  the  lossless  case,  where  impulse  functions  are  included  in  the  transform  of  Eq.  (3.7). 

The  efficient  evaluation  of  the  integral  in  Eq.  (3.2)  is  quite  difficult  and  of  prime  importance. 
The  inner  product  integral  is  of  the  form 


dkz 


(3.9) 


where  i,  j  €  xj  and  F(.kz)  and  Y,j  have  no  real  poles.  The  two  integration  contours  used  are 
shown  in  Fig.  33.  Figure  33a  is  for  positive  traveling  waves  and  Fig.  3.3b  for  negative  traveling 
waves.  The  pole  is  located  at  the  propagation  constant  for  a  given  region.  A  good  discussion  of  the 
contour  selection  for  many  mode  problems  is  given  in  Collin  and  Zucker  [20].  There  are  an  infinite 
number  along  the  imaginary  axis  due  to  the  evanescent  modes.  The  location  of  these  cannot  be  deter¬ 
mined  simply.  Therefore,  the  integration  contour  cannot  be  closed.  It  is  then  necessary  to  evaluate 
the  principal  value  of  the  integral  numerically.  Davis  and  Rabinowitz  [21]  discuss  methods  of 
approximate  integration  over  infinite  intervals.  Unfortunately,  the  techniques  generally  used  for 
such  integrals  are  not  readily  applied  to  this  case.  The  kernel  has  a  two-dimensional  variation  which 
makes  the  application  of  asymptotic  techniques  very  difficult.  A  satisfactory  truncation  point  for  the 
integral  must  be  selected  to  give  an  acceptable  error.  The  method  of  foldover  [22]  may  be  applied. 
This  technique  may  be  summarized  by 


J  Fi(kz)dk:  =  f  F  ,(*.  )dk.  +  f  F^k.'kik.  +  j  w  £es(-fi  ,X 

— OO  -“OO  —  fi  j 

-W,  0 

=  /  F,(*.)d*z  +  /  ^(-20,-*,) +  F  ,(*,)</** 

-0i  1 

oo 

+  f  F  {(kz  )dkz  +  j-rr  Res(—0t). 


(3.10) 


oo  oo 

f  Ffkz)dkz  =  f  Ffkz)dk.  +  /  F  fkz)dkz  -  jir  Resi-fiJ 

-OO  -OO  /}j 

0  -*«  [  . 

=  f  F  fkz)dkz  +  /  [Fi(Jkt)  +  Fj(231-*i)|j*2 


i7  /A:z  )d*2  —  jir  Iks  (— 3 ,). 


(3.11) 


The  new  kernel  becomes  bounded  at  the  location  of  the  singularities  in  the  old  kernel.  This 
allows  a  more  accurate  approximation  near  the  singularity  and  an  increase  in  the  speed  of  the  compu¬ 
tation.  The  kernel  is  oscillatory  in  nature,  which  increases  the  difficulty  of  the  numerical  integration. 


33.  Numerical  Results  for  the  Single  Discontinuity 

Ex  is  the  dominant  field  component  in  the  slot.  It  is,  therefore,  a  satisfactory  approximation  to 
set  Ez  =  0.  Equation  (3.6)  then  becomes 


<i»f~Ex>  =  0 


(3.12) 


for  £m  orthogonal  to  J  x . 


Now  consider  a  single  pulse  approximation  for  Ex  (x  ),  90  that  Ex  can  be  written  as 


(3.13) 


Ex  +  ateJI>'z)P,(s,wl)(l-u_l(z)) 

+  a2e~J$*  PICs,w2)u_,(r) 

where  fij  and  /32  are  the  propagation  constants  for  the  dominant  mode  for  z  <  0  and  z  >  0,  respec¬ 
tively.  The  testing  functions  used  are 

£ ,(x  ^  )  =  Px  is  ,w ,)  Pz  ( e  J*'z 

*  Fx(s,wa)P,(-^pa)«',*a*.  (3.14) 

The  integral  in  Eq.  (32)  has  been  evaluated  numerically  using  both  the  IMSL  routine  DCADRE  and 
the  trapezoidal  scheme.  A  more  efficient  computer  program  can  be  written  using  the  trapezoidal  rule. 
Integral  limits  of  ±  15  with  nineteen  spectral  terms  in  the  Jkx  summation  prove  to  be  adequate  for 
convergence.  A  trapezoidal  integration  with  700  points  is  suitable.  The  foldover  method  is  used  in 
the  integral  evaluation. 

The  scattering  parameters  of  the  Junction  computed  by  this  procedure  are  not  stable  with  respect 
to  variations  in  the  pulse  length  a.  While  somewhat  reasonable  results  can  be  achieved  for  certain 
cases,  the  approach  of  testing  a  semi-infinite  unknown  with  a  finite  pulse  does  not  give  reliable 
results. 

If  loss  is  introduced  into  the  basis  functions  of  Eq.  (3.13)  and  Galer kin’s  method  is  used  in  Eq. 
(3.12),  the  poles  move  off  the  real  kt  axis.  The  unknowns  are  tested  over  the  entire  semi-infinite 
region.  The  scattering  parameters  vary  significantly  as  a  function  of  the  loss  coefficient. 

3.4.  Theoretical  Analysis  for  an  Infinite  Periodic  Array  of  Discontinuities 

A  unit  cell  in  an  infinite  array  is  shown  in  Fig.  3.4.  The  spectral  relationship  of  Eq.  (3.1)  can  be 
applied  to  this  geometry.  However,  now  kt  takes  on  discrete  values. 
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A  discussion  of  periodic  structures  is  given  in  Collin  and  Zucker  [20].  Floquet’s  Theorem  may 
be  used  to  relate  the  field  quantities  at  z  =  —a  and  z  -  a.  The  fields  may  be  expanded  in  a  Fourier 
series  in  the  z-direction.  For  the  E-field, 


Ef(x,yj)  =  Ef„(x,y)e~J0nZ 
1  “12 

Eqnix,y)=—  f  Ell(x,yrz)ejfi*zdz 

a  -or/2 


where  q  €  \xj) 


fin 


=  fio 


2n  it 
a 


(3.15) 


In  the  application  of  Eq.  (3.1)  to  this  geometry,  field  quantities  in  the  plane  of  the  fin 
(y  =  1 1  +  2d)  are  considered.  Ex  and  Ez  in  this  plane  can  be  expanded  in  terms  of  appropriate  func¬ 
tions,  and  a  homogeneous  equation  of  the  form  Ax  -  0  can  then  be  obtained  by  using  Galerkin’s 
method. 

Consider  a  representation  for  Ex  and  E.  given  by 


Ex(xj)  =  $.(*.*) 

i=i 

/ 

Ez(xj)=  5}  by  T)j(xj).  (3.16) 

The  weighting  for  the  discrete  spectral  terms  is  given  by 


a/2  2b 

\(kx3„)  =  •—  f  f  ;i(xa)ejtl*x*li'2,dxdz 

a  —or/2  0 

ar/2  2b 

f)Mxfin)- -rr—  f  f  t )j(xj)eMxX+*HZ)dxdz.  (3.17) 

®  —a/2  0 


The  homogeneous  equation  becomes 
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E  <lP'  E >  +  E  <1 P .  y*z  £  M,  >  =  °-  p  =  w 

n  t =1  n  J=l 

E  E*,L>  +  E  <*>,.!%*  EM;>  =°-  <?  =  w  (3.i8) 

it  i=l  n  j  — 1 

The  solution  of  det(A)  -  0,  where  A  is  the  operator  matrix  in  Eq.  (3.18),  gives  the  solution  for  0„ 

Consider  the  Gfc0  a ,  0a)  diagram  shown  in  Fig.  3.5.  Several  coupled  modes  are  shown  in  this 
diagram.  Stopbands  occur  when  0a  equals  an  integral  multiple  of  tr.  When  0a  =  ir,  coupling  occurs 
between  the  n  *  0  Floquet  mode  and  the  n  *  *1  mode.  That  is,  there  is  coupling  between  a  forward 
traveling  wave  and  a  negative  traveling  wave.  Note  that  the  example  in  Fig.  3.5  is  for  a  TEM  wave 
in  a  periodically  loaded  media  which  has  no  low-frequency  cutoff. 

The  analysis  outlined  here  is  implemented  numerically  in  the  next  section.  This  gives  an 
approximation  to  the  ( ke  a ,  0a)  curve  over  a  particular  range.  The  range  of  interest  here  is 

0  <  0a  <  tt. 


3.5.  Numerical  Results  for  an  Infinite  Periodic  Array  of  Discontinuities 

The  0to  a,0a)  diagram  for  an  infinite  periodic  array  of  discontinuities,  as  shown  in  Fig.  3.1(b),  is 
determined  in  this  section.  The  basis  functions  used  to  approximate  the  fields  in  Eq.  (3.16)  are 

£i(jc,z)  =  Px  is  ,w  j)  Pz  (0/») 


&*(x,z)  =  Px(s,w ,)  Pt (0,a)  —  Pt (Of, ) J  sgn (x) 


t)|(x  £ )  =  2~ 


,  (x  —  s ) 


w 


(3.19) 


The  computed  results  for  a  portion  of  the  (*<>a  ,  0a)  diagram  are  shown  in  Fig.  3.6  together  with  the 
fin-i in*  dimensions  for  this  example.  The  approximation  used  for  the  slot  fields  results  in  some  spun- 
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ous  solutions  near  fid  =  it  and  poor  sensitivity  to  discontinuity  variations. 

3£.  Concluding  Remarks 

This  chapter  used  a  two-dimensional  spectral  relationship  in  the  fin  plane,  between  the  slot 
fields  and  the  fin  currents.  The  moment  method  was  then  used  to  solve  for  the  slot  fields.  Both  the 
problem  of  a  single  step  discontinuity  and  an  infinite  periodic  array  of  discontinuities  were  studied. 
In  the  case  of  the  single  discontinuity,  efficient  evaluation  of  the  inner  product  integrals  is  important 
to  minimize  computation  time. 


CHAPTER  4:  ITERATIVE  AND  VARIATIONAL  SOLUTIONS  FOR  FIN-LINE 


DISCONTINUITIES 


4.1.  Introduction 

This  chapter  considers  various  iterative  solutions  for  the  fin-line  discontinuity  problem.  The 
discontinuity  considered  in  this  chapter  is  that  of  Section  32  which  is  shown  again  in  Fig.  4.1.  The 
emphasis  of  this  chapter  will  be  on  two  iterative  procedures  which  are  applied  to  a  suitable  formula¬ 
tion. 

Solution  of  electromagnetic  problems  using  iterative  techniques  is  advantageous  when  there  are 
a  large  number  of  unknowns.  In  fact,  when  there  are  a  very  large  number  of  unknowns,  it  is  the 
only  way  to  find  a  solution.  If  the  number  of  unknowns  is  such  that  the  matrix  can  be  inverted  on 
the  computer  being  used,  the  moment  method  may  be  preferable. 

The  discontinuity  problem  is  formulated  using  unknowns  in  the  junction  plane,  in  Section  42. 
The  solution  procedure  using  a  conjugate  gradient  iterative  scheme  is  outlined  in  Section  43.  Numeri¬ 
cal  results  for  the  conjugate  gradient  method  are  given  in  Section  4A  Solution  for  the  junction  plane 
unknowns  using  the  generalized  variational  procedure  is  outlined  in  Section  4.5,  and  numerical 
results  are  given  in  Section  4.6.  Conclusions  are  drawn  in  Section  4.7. 

42.  Formulation  using  Unknowns  in  the  Junction  Plane 

A  waveguide  junction  problem  can  be  formulated  using  the  mode  functions  either  side  of  the 
discontinuity  [23].  Consider  the  general  junction  of  Fig.  42.  The  transverse  field  can  be  expressed  in 
terms  of  the  mode  functions  in  guides  A  and  B.  Consider  I  modes 
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ea  e^**  +  E  a, e«  e***  2  <  0 

—  1*1 

E,  =  ;  . 

E  6,  **  *">4,z  *  >  0 

i»l 


(4.1a) 


hm  e^*1  ~  E  OiK  ey*z  2  <  0 

_  i*  1 

=  /  . 

T.bihtt  e^1  2  >  0 


(4.1b) 


Rather  than  pursuing  the  regular  mode-matching  formulation,  consider  an  alternate  approach  as 
follows.  Figure  4J  shows  an  equivalent  representation  for  the  junction,  employing  equivalent  mag¬ 
netic  currents  over  the  apertures.  The  total  transverse  field  for  the  equivalent  problem  is 


a,  = 


em  y<*2  —  e 

*-z  +  E4 «v-z 

«-l 

2  <0 

i*  1 

2  >  0 

hm  e^z  +  hm 

ey-z  -t*i  K  «y*z 

2  <  0 

»«1 

2  >  0 

(4.2a) 


(4.2b) 


The  fields  produced  in  guide  A  by  K  are  Ea  {K ),  Ha  (K ),  and  in  guide  B  by  —K  are 
Eb  ( —K  ),  Hb  {—K  ).  For  continuity  of  E, , 


K  *  2XEa<lK  )l2^j  =  E  di  zxe^ 

i*l 

K  ■  rx£j(— iT )ljm  =  J  i|  2 Xfj, . 

i»I 

Continuity  of  H,  requires 
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2/Ta,  =  jld.K  +  'tb.K.  (4.4) 

i=i  (=i 

The  solution  for  the  unknown  quantity,  K ,  could  be  obtained  via  several  numerical  approaches. 
The  solution  for  k  via  certain  iterative  techniques  is  of  interest  in  this  chapter. 

Taking  inner  products  of  Eq.  (4.3),  and  using  the  orthogonality  of  Eq.  (2.18)  with  normalized 
mode  functions. 


i 

E 

i  =1 

=  dt 

I 

E 

i»i  s 

=  6*.  (4.5) 

Substitution  for  the  d’s  and  b’s  in  Eq.  (4.4)  results  in  an  integral  equation  which  may  be  used  to 
solve  for  the  unknown  operative  magnetic  current  K . 


ff  K  ds  =  E4  ff  '  2  Xea  ds 


f f  hht  ■  K  ds  -  E*>,  ff  h^ 


2  (i  ,y )  *  t  f  f  hm(*,?)-  K(*,y)d*dy 

«*1  s 

+  E  K  (*  .y )  //  ft*  (x-  ,y )  •  k  ,y  My.  (4.6) 

i*l  s 

This  equation  may  be  interpreted  in  terms  of  the  Green’s  function  for  each  guide  as  follows. 
The  dyadic  Green’s  function  can  be  found  simply  from  the  mode  functions  using  a  procedure  outlined 
in  Collin  and  Zucker  [241  Equation  (4.6)  can  be  written  as 


47 


2 h,«  //  [C„  +  G*]UM<fc' 


where 


h  — 

(X,y) 

n,  - 

hmy  (x,y) 

X*  (*.y) 

K  = 

Ky  (x,y) 

E  ^»«i  (*  .y  )*,«  (*\/)  E  V*  k  «y  (*  .y1 ) 

<  =*1  i*l 

E  h<*j  C*  .y  )*7o  (*\y )  E  ^7  (*  .y  (*  ,y ) 

»=1  i*l 

q  €  {a>}. 


(4.7) 


For  the  fin-line  discontinuity  of  Fig.  4.1,  the  aperture  is  the  cross-section  of  the  shielding  enclo¬ 
sure.  The  number  of  modes,  I,  in  each  guide  will  be  small  due  to  the  difficulty  of  finding  a  large 
number  of  orthogonal  modes. 

Following  the  solution  for  K ,  the  mode  coefficients  can  be  found  using  Eq.  (4.5).  The  scattering 
parameters  can  then  be  determined  from  these  mode  coefficients  following  a  re-normalization  for  unit 
power.  Assuming  that  the  modes  have  been  normalized  for  f  f  en  Xhn  ds  =  1,  the  scattering  param- 

s 

eters  for  the  single  propagating  mode  are  given  by 


Sn »  a  II  N  i  ‘  K  ds 


f  f  *«*Kds 


1/2 


1/2 


q  €  {ajt>} 
n  €  {1,2}. 


Iff  e^Xh^ds 


(4.8) 
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43.  Conjugate  Gradient  Solution  for  Unknowns  in  the  Junction  Plane 

The  conjugate  gradient  method  [25H28]  has  been  in  existence  for  some  time.  However,  it  is  only 
recently  that  it  has  been  applied  to  the  solution  of  electromagnetic  problems  [27].  The  method  was 
developed  by  Hestenes  and  Stiefel  [25]. 

For  n  unknowns,  the  method  will  converge  in  n  steps,  if  no  rounding-off  errors  occur.  The 
method  is  quite  simple.  At  each  step,  the  estimate  is  an  improvement  over  the  previous  step,  and  the 
process  can  be  restarted  at  any  step  using  the  last  value  as  the  new  initial  estimate. 

Consider  the  linear  equation  L  f  =  g,  where  L  is  a  general  non-singular  matrix.  LAL  is  a 
positive  definite  operator,  but  is  in  general  non-symmetric.  The  system  L  f  =  g  is  equivalent  to 

LALf  =  LAg.  (4.9) 

Another  equivalent  equation,  which  results  in  a  symmetric,  but  not  a  positive  definite  operator,  is 
LT Lf  =  LT g.  The  conjugate  gradient  technique  may  be  applied  to  Eq.  (4.9).  The  iterative  scheme 
is 

ra=  g  ~  Lf0 

Po  -  LAr0 

=  IILAri112 
*  IILp.  ll2 

/i+i  =  /.  +«iPi 

'•,+1  =  n  -aiLpi 

b  IILA>W|2 

i_  IILAr,  II2 

A+i  =  LAri+I  +biPi. 

The  inner  products  in  this  scheme  are  defined  such  that 
</X/  >  =  <LA  f ,/  >. 


(4.10) 


/  II 2  =  <  /  ,/  >  and 


An  estimate  / 0  of  the  unknown  /  is  used  to  generate  the  residual  r0  and  the  direction  p0 .  In 
each  cycle,  following  the  determination  of  the  estimate  /,,  the  residual  r,  and  the  direction  p,,  the 
next  estimates  for  /1+,/-,+1,pt+j  are  computed.  The  residuals  r0  r ,,  ■  •  are  mutually  orthogonal,  and 
the  direction  vectors  potp „ ••  •  are  mutually  conjugate,  therefore 

<r,Vy  >  —  0  l*j  (4.11) 

<PiJLpj>  =0. 

Each  estimate  /i+1  will  be  closer  to  the  solution  than  The  norm-squared  of  the  residual  can  be 
used  as  a  measure  of  the  error 

erri  —  II  Tj  II2.  (4.12) 

However,  this  is  not  alway  reliable,  as  II  II 2  may  not  decrease  at  each  step. 

Consider  now  the  fin-line  discontinuity  problem  of  interest,  which  is  shown  in  Fig.  4.1.  This 
problem  was  formulated  in  Section  42.  The  conjugate  gradient  method  can  be  applied  to  Eq.  (4.6)  to 
solve  for  K  (x  ,y ).  The  unknowns  are  Kx  and  Ky  over  a  grid  of  points  in  the  (x,y)  plane  at  z  -  0. 
An  (mj i)  set  of  grid  points  represents  2 mn  unknowns.  A  reasonable  approximation  for  the  unk¬ 
nown  requires  a  large  number  of  unknowns. 

The  inner  product  for  this  problem  is  given  by 

<x«y  >=//*>  d*dy-  (4.13) 

S 

4.4.  Numerical  Results  for  the  Conjugate  Gradient  Solution  in  the  Junction  Plane 

The  mode  functions  are  found  using  the  procedure  outlined  in  Chapter  2,  resulting  in 
en  (kx  ,y)  ,h„  (*x  ,y  )•  For  convenience,  many  of  the  processes  in  the  iterative  routine  of-  Eq.  (4.10)  are 
carried  out  in  the  spectral  domain.  The  Fast  Fourier  Transform  is  utilized  throughout  for  transfor¬ 
mations  in  the  x-direction.  Although  there  are  2 mn  total  unknowns  (for  an  (m  ji )  array  of  points). 


50 


(mji)  arrays  may  be  used  throughout  the  computer  program. 

As  indicated  in  Chapter  2,  it  is  difficult  to  obtain  a  large  set  of  orthogonal  mode  functions  for 
use  in  the  Green’s  function.  This  necessitates  computation  with  as  few  as  three  modes. 

The  limiting  case  of  a  uniform  slot  is  a  good  example  to  consider  first.  A  relative  error  will  be 
used  as  a  figure  of  merit,  where  the  residual  is  normalized  to  the  incident  field. 


reL  err,  — 


(4.14) 


Figure  4.4  shows  the  results  at  each  iteration  using  a  poor  initial  estimate  of  kx  =  ky  -  0  +  j  0. 
Notice  the  rapid  convergence  to  the  solution.  This  data  is  for  m  =17,  n  =33  w,  =w2  =  1mm  and 
three  modes.  There  are  therefore  1122  unknowns.  Using  a  DEC  VAX  11/780  computer,  ten  itera¬ 
tions  using  this  array  size  require  approximately  sixteen  seconds  cpu  time.  This  time  can  be  reduced 
somewhat  by  using  a  radix-two  FFT  and  zero  padding. 

Applying  the  same  procedure  to  the  case  of  differing  slot  widths  on  either  side  of  the  junction 
does  not  yield  such  satisfactory  results.  Figure  4.5  shows  the  results  for  m  =17,  n  =33, 
w ,  =  1mm,  w  2  —  2mm  and  three  modes,  s  u  goes  to  one  because  the  iterative  scheme  finds  a  solution 
which  corresponds  to  a  higher-order  mode,  and  is  orthogonal  to  the  Green’s  function  modes  in  one 
region.  The  conclusion  is  that  with  the  small  number  of  modes  which  can  realistically  be  found  and 
the  large  number  of  degrees  of  freedom  in  the  unknown,  this  approach  will  not  yield  satisfactory 
results.  It  is  necessary  to  have  the  number  of  unknowns  similar  to  the  number  of  modes  in  the 
Green’s  function.  Similar  results  are  obtained  when  the  number  of  grid  points  is  quadrupled. 


4. 5.  Solution  in  the  Junction  Plane  Using  the  Generalized  Variational  Technique 

The  approach  presented  here  also  considers  unknowns  in  the  (x,y)  plane  of  the  junction.  Again 
consider  Eq.  (4.6),  which  is  of  the  form  Lf  -  g ,  where  L  is  an  integral  operator,  /  is  the  unknown, 
and  g  is  the  known  incident  field. 


I  tarot  Ion  number 


Relative  error  for  the  conjugate  gradient  solution  of  a  uniform  fin-line  with  a  WR-28 
shield.  1 1  -  1 2  "  3.429  mm,  2d  -  0.254  mm,  er  -  2 22,  2b  -  2.556  mm,  w ,  -  w2  -  1 
Tnm,  f  -  30  GHz,  m  -  17,  n  -  33.  The  resultant  scattering  parameters  are:  S  u  -  1.192 
x  10"5  +  j  1.0486  x  10“*,  S  a  -  1.000  +  j  1.049  x  10"*. 


I  t  «r a  t I  on  number 


Relative  error  for  the  conjugate  gradient  solution  of  a  discontinuity  in  fin-line  with  a 
WR-28  shield-  Z,  -  l2  -  3.429  mm,  2d  -  0.254  mm,  er  -  2.22,  2b  -  2.556  mm, 
w ,  =  1  mm ,  w  2  -  2  mm,  /  *  30  GHz,  m  -  17,  n  ■  33.  The  resultant  scattering  param¬ 
eters  are:  Sn  -  1.000  +  j  6.487  x  10-3  S2l  -  -1.904  x  10-3  +  j  1.636  x  10-2 


The  unknown  quantity  x  may  be  expressed  as  a  superposition  of  characteristic  functions  which 


are  to  be  determined.  /  can  be  expressed  as 


/  =  E  cnfn. 


n  *0 


(4.15) 


The  /„  are  L-orthogonal,  and  satisfy 


<fn^fm>-T)nKm  (4.16) 

where  is  the  Kronecker  delta  function.  A  suitable  set  of  characteristic  functions  can  be  generated 

by 


fn 


II “n  H2 


+  fn- 1 


where  u„  is  an  auxiliary  function  which  satisfies  the  orthogonality  property 


(4.17) 


^  ^  in  • 

The  auxiliary  function  can  be  found  by 


(4.18) 


=  a*-,  ~  $„-i  -Lx„-i  (4.19) 

where 

Vn-l 

From  Eqs.  (4.15)  and  (4.16), 

Cn  ~  ~  ./„•*>•  (4.20) 


The  generation  routine  can  be  summarized  as 
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fo  = 


■Hn  =  <fn>Lfn  > 

.  _  </„*> 


“n+l  =“n  "  ~  L/" 

’ll. 

f  =  fflll  4.  / 

<  «+l  „  mi  ^  J  n. 

IlKn+lll2 


(4.21) 


The  generalized  variational  technique  may  also  be  applied  to  the  equation  LALf  -  LA g, 
where  now  the  generation  routine  becomes 


Vn  =  <fnfLALfn> 

<fn*Ag> 

c" - z - 

’ll 

Un+l  =  u„  -±-LALfr 

f  m  f 
J  n  + 1  „  n  i  '  J  n. 


(4.22) 


These  routines  can  be  used  to  solve  for  the  unknown  magnetic  current  in  the  junction  plane  of 
the  fin-line  discontinuity.  A  good  choice  for  /„  is  zxea  „  where  ea ,  is  the  incident  transverse  E-field 

A 

mode  function,  with  a  possible  alternative  being  z  x  (  ea ,  +  eb  ,  ).  The  coefficients  c„ ,  and  successive 
characteristic  functions  /„  are  generated  by  Eq.  (4.21)  or  Eq.  (4.22). 


4j6.  Numerical  Results  for  the  Generalized  Variational  Solution  of  the  Fin-Line  Discontinuity 
Problem 


The  results  in  this  section  are  generated  using  three  modes  in  the  Green’s  function.  Two  basis 
functions  for  Ex  and  Et  are  used  to  generate  the  modes.  Referring  to  Eqs.  (2.8)  and  (2.9),  these  func¬ 
tions  are  for  p  «  0, 1  and  q  *  2, 4.  It  is  important  to  verify  that  the  numerical  results  satisfy  the 
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following  relationship  for  the  single  propagating  mode. 

IS„I2  +  IS2,I2  =  1.  4.23 

Referring  to  Eq.  (4.21)  and  Eq.  (4.22),  u0  is  chosen  as  2Xea ,.  Two  different  unilateral  fin-line 
geometries  are  considered.  There  is  a  shortage  of  data  from  other  authors  for  comparative  purposes. 

Figures  4.6-4.9  show  the  computed  scattering  parameters  for  a  step  discontinuity  in  fin-line 
with  a  WR-28  shield.  The  equation  Lf  —  g  is  used  with  one  and  two  characteristic  functions.  A 
monotonically  convergent  solution  as  a  function  of  the  number  of  characteristic  functions  is  not 
guaranteed.  A  comparison  with  the  mode-matching  data  in  [17]  is  given.  However,  it  appears  that 
Eq.  (4.23)  is  not  satisfied  well  in  [17],  so  definite  conclusions  cannot  be  drawn  regarding  the  accuracy. 
A  comparison  between  the  equations  Lf  =  g  with  one  and  two  characteristic  functions  and  the 
equation  LALf  —  LAg  with  two  characteristic  functions  is  given  in  Figs.  4.10-4.13. 

The  scattering  parameters  for  a  step  discontinuity  in  fin-line  with  a  WR-62  shield  are  shown  in 
Figs.  4.14-4.17.  The  equation  Lf  —  g  is  used  here.  Equivalent  circuits  for  this  problem  using  S  u 
have  been  generated  by  ElHennawy  and  Schunemann  [16]  using  a  mode  matching  technique.  A  com¬ 
parison  with  these  results  is  given.  An  indication  of  the  modal  solution  accuracy  and  the  satisfaction 
of  Eq.  (4.23)  cannot  be  determined  from  the  work  of  ElHennawy  and  Schunemann  [4i  [16].  An 
approximate  convergence  criterion  for  the  maximum  number  of  characteristic  functions  p  with  q 
modes  in  the  Green’s  function  q  is  p  <  q  —1. 

47.  Concluding  Remarks 

This  chapter  considered  two  procedures  for  evaluating  unknowns  in  the  cross-sectional  plane  of 
the  fin-line  discontinuity.  The  scattering  parameters,  or  equivalent  circuit  for  the  discontinuity,  may 
be  found  from  these  unknowns. 
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The  conjugate  gradient  method  was  considered  first.  It  did  not  give  satisfactory  results  for  the 
discontinuity  problem.  This  was  because  the  number  of  degrees  of  freedom  in  the  unknown  was 


*x  % 


magn i tud*  a  1 t 


phas*  «2t  CdttO- 


•Fr««ju«ncy  Chz? 


Figure  4.9. 


PhaaeCS  21)  for  a  discontinuity  in  fin-line  with  a  WR-28  shield.  l\mlt  “.3.429  mm,  2d 
-  0.254  mm,  er  -  222,  2b  -  2.556  mm,  w ,  -  1  mm,  w  2  -  2  mm,  m  -  17,  n  -  33. 
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Figure  4.10. 


IS  u  1  for  a  discontinuity  in  fin-line  with  a  WR-28  shield.  I ,  «  1 2  m  3-429  mm,  2d 
0254  mm,  €r  -  22.2,  2 b  -  2.556  mm,  w,  -  1  mm,  w2  -  2  mm,  m  -  17,  n  -  33. 
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Figure  4.11. 
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Phase  (S  u)  for  a  discontinuity  in  fin-line  with  a  WR-28  shield.  I ,  - 1 2  •  3.429mm,  2d 
-  0.254  mm,  6r  «  2.22,  2b  “  2^556  mm,  w,»l  mm,  w  2  -  2  mm,  m  -  17,  n  -  33. 
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Figure  4.13.  Phased  21)  for  a  discontinuity  in  fin-line  with  a  WR-28  shield.  3.429  mm,  2d 

”  0.254  mm,  er  m  2.22,  2b  -  2.556  mm,  w  j  -  1  mm,  w  2  -  2  mm,  m  -  17,  n  -  33. 
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Figure  4.14.  IS  u  I  for  a  discontinuity  in  fin-line  with  a  WR-62  shield.  1 1  -  l2  m  7.773  mm,  2d 
0254  mm,  €.  -  22,  2b  -  3.95  mm,  w ,  -  3  mm,  m  -  17,  m  -  33. 


Figure  4.15.  Phase  (5  u)  for  a  discontinuity  in  fin-line  with  a  WR-62  shield.  I ,  ■  l_2  -  7.773  mm, 
2d  -  0254  mm,  €r  -  22,  2b  -  3.95  mm,  w ,  -  3  mm,  m  -  17,  n  -  33. 
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much  greater  than  terms,  or  modes,  in  the  Green’s  function.  It  is  not  realistic  to  find  a  large  number 
of  modes  for  the  Green’s  function.  The  uniform  guide  solution  was  reconstructed  very  efficiently. 

The  other  approach  outlined  was  the  generalized  variational  technique,  which  involved 
representing  the  junction  plane  magnetic  current  in  terms  of  characteristic  functions.  The  routine  for 
generating  these  functions  and  determining  their  coefficients  was  outlined.  The  results  using  this 
approach  were  satisfactory. 

The  difficulty  of  finding  an  accurate  large  set  of  mode  functions  is  a  limitation  of  the  approaches 
covered  in  this  chapter.  For  this  reason,  it  is  advantageous  to  look  at  alternate  formulations  which  do 
not  require  the  generation  of  these  mode  functions. 

Formulating  the  problem  in  the  fin-plane  avoids  the  need  for  computing  the  mode  functions. 
The  Green’s  function  determined  in  Chapter  2  can  be  used,  as  it  was  in  Chapter  3  for  the  moment 
method  analysis.  The  author  has  investigated  several  iteration  procedures  which  may  be  used  in  con¬ 
junction  with  such  a  formulation.  This  work  involved  a  variant  of  the  Spectral  Iterative  Technique 
(SIT)  [291  For  guaranteed  convergence,  it  would  be  advantageous  to  apply  the  conjugate  gradient 
technique  to  this  problem.  These  approaches  are  worthy  of  further  consideration. 


CHAPTER  5:  DIELECTRIC  WAVEGUIDE  FILTERS 


5.1.  Introduction 

Dielectric  waveguide  (DWG)  filters  are  of  particular  interest  for  millimeter- wave  end  optical 
applications.  A  number  of  applications  have  been  reported  recently  in  the  millimeter-wave  and  quasi- 
optical  areas  [30H321  and  optical  grating  filters  have  been  reported  in  [33].  Grating-type  dielectric 
waveguide  filters  have  several  advantages  over  alternative  configurations  such  as  the  ring-resonator 
filter  [34],  [35].  In  particular,  the  ring  resonator  should  be  a  number  of  wavelengths  in  circumference 
for  satisfactory  performance,  which  implies  closely  spaced  spurious  pass-  or  stop-bands.  The  filter  can 
easily  be  incorporated  into  an  integrated  system  and  may  be  realized  by  a  series  of  discontinuities  such 
as  surface  or  dielectric  variations. 

There  has  been  some  work  done  on  mode-matching  methods  [36]  in  an  effort  to  analyze  discon¬ 
tinuities  in  a  two-dimensional  dielectric  waveguide.  It  is  difficult  to  apply  these  techniques  to  the 
analysis  of  grating  filters  considered  in  this  chapter.  It  is  proposed  that  a  simple  transmission  line  may 
be  used  to  analyze  such  a  grating  structure.  Experimentation  has  verified  this  approach.  The  theory 
developed  is  valid  for  a  number  of  dielectric  waveguide  structures. 

Section  3.2  gives  details  of  the  theoretical  analysis  of  the  dielectric  waveguide  grating  filter. 
Experimental  results  are  given  in  Section  5-3.  Conclusions  are  outlined  in  Section  5.4. 

52.  Analysis  of  Dielectric  Waveguide  Crating  Structure 

A  stepped  dielectric  grating  structure  and  transmission  line  model  are  shown  in  Fig.  5.1.  The 
waveguide  known  as  image  guide  (rectangular  dielectric  guide  on  a  ground  plane)  will  be  considered 
specifically,  although  the  approach  may  be  generalized  to  other  dielectric  waveguides.  It  is  proposed 
that  this  model  may  be  used  for  an  approximate  analysis  of  the  grating  filter. 

The  dielectric  waveguide  supports  hybrid  modes  which  may  be  expressed  as  a  sum  of  longitudinal 
section  electric  (LSE)  and  longitudinal  section  magnetic  (LSM)  modes.  The  LSE  fields  may  be  expressed 


vv> 


(5.1a) 


E  =  -  j  o>  V  X  IX 

fi=«^n»+vv'a  (5.1b) 


and  the  LSM  fields  as 


E  =  *<?IL  +  V(er*  V  •  IL)  (52a) 

H  =  j  at  «oV  X  IX  (52b) 


where  IX  and  JX  are  the  magnetic  and  electric  Hertzian  potentials,  respectively. 

With  IX,  =  y  <(>h(x,y)e~jet  =  y0A  and  IX  =  y4>,  the  dielectric  waveguide  fields  may  be 
expressed  in  terms  of  electric  and  magnetic  scalar  potential  functions: 
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The  effective  dielectric  constant  method  [37l  [38]  is  a  very  suitable  approach  for  finding  the  propagation 
constant  in  planar  dielectric  waveguides.  Various  forms  of  dielectric  waveguides  may  be  analyzed  in 
this  manner.  Consider  the  image  guide  structure  of  Fig.  52.  Using  the  effective  dielectric  constant 
method  and  matching  the  fields  in  each  region,  the  following  eigenvalue  equations  may  be  found.  The 
equation  for  k7  is 

tr  V*<?(«r  ~  1 )  —  V  cos  k,y ,  —  sin  k7y ,  =  0  (5.4) 


where  er  is  the  relative  permittivity  of  the  guide.  After  solving  for  k.,  the  effective  dielectric  con- 
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Figure  52. 


w  =  x2-x1 


Image  guide  with  coordinate  system. 
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stant  (ee2)  for  region  2  may  be  found: 


The  equation  for  kx  is 

-  1)  —  2*123sin  (x2  -  x ,)  + 


+  2kx  >/k$(et 2  -  l)  -  kx2  cos  kx  (x2  -  x  „  =  0. 

Thus  the  guide  propagation  constant  becomes 

0  =  V«.2*<?  ~kx2.  (5.7) 

The  dominant  mode  will  be  considered  in  the  following  analysis.  The  waveguide  characteristic 
impedance  may  be  defined  as 

Z  =  Vo  J-  (5.8) 

where  tjo  *0  are  the  f ree-space  wave  impedance  and  propagation  constant,  respectively. 

In  order  to  find  the  filter  transfer  function,  consider  the  unit  cell  shown  in  Fig.  53.  The  ABCD 
matrix  for  the  unit  cell  may  be  obtained  by  multiplying  the  matrices  for  two  line  sections  [39].  The 
unit  cell  matrix  is 


A  —  co^0c Iq ) ccrtBuc lffC )  —  ■  -qG  sm(j8c  Iq )sm(0jvC lNG ) 

Pc 


’  O  _ _ /  a  »  ^  0  .  _ 


B  —  jt)o  „  ccrtPc Zg ) ®n(0G Zg )  +  "3 —  sin(0G Z g ) cos (0yc ^nc ) 

PVG  PC 


C  *  ~r~  sin(0G  Zg  )  cos  (0c  Z#g  )  +  cos  (0o  Iq  )  sin  (0c  (vc  ) 

7)o  *o  *o 


^  *  cos(0g  Zc )  c*  (0w  Ing  )  ”  sin(0G  Zg)  sin  (0tfG  Zjvg  )• 


The  ABCD  matrix  for  a  grating  structure  of  n  unit  cells  is  obtained  by  raising  the  matrix  for  a  single 
cell  to  the  nth  power. 

The  transducer  loss  ratio  is  given  by  [39] 


4 ScRl 
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(5.10) 


where  Pi  is  the  power  delivered  to  Rl  and  =  I  Vg  I  2/4Rq  is  the  available  power  from  the  genera¬ 
tor.  This  results  in  a  method  for  computing  the  filter  insertion  loss  (Pi  !Pmail')  as  a  function  of  fre¬ 
quency.  Although  this  technique  is  approximate,  the  response  predicted  from  such  an  analysis  agrees 
well  with  measurements. 

Fields  in  a  periodic  structure  can  be  expressed  in  terms  of  spatial  harmonics,  according  to  Floquet’s 
Theorem  [l9l  The  propagation  constants  are 


Pn  =  0o  +  2~~,  n  =0,±1,±2. 


(5.11) 


where  d  is  the  grating  period.  It  is  useful  to  look  at  the  dispersion  curves  on  a  k&IJ}  d-diagram,  such 
as  Fig.  3.6.  Coupled  mode  theory  may  be  used  to  explain  the  filter  characteristics.  There  is  a  stop  band 
when 


0  d  =  it 


(5.12) 


where  0  is  the  propagation  constant  in  the  grating.  At  this  frequency  there  is  coupling  between  the 
spatial  harmonics  (0q)+  and  (0oi)—  which  propagate  in  opposite  directions.  The  waves  are  guided  by 
the  structure  (slow-wave  region).  The  grating  exhibits  a  reflection  coefficient  close  to  one  in  the  stop- 
band,  and  close  to  zero  on  either  side. 

Consider  ports  n  and  n+1  in  the  unit  cell  of  Fig.  5-3,  when  the  cell  is  part  of  an  infinite  periodic 
structure. 
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(5.13) 


For  a  lossless  reciprocal  network  [40], 


«*  3  d  =  (5.14) 

The  following  relationship  may  thus  be  derived 


cos  fid  =  cos  &clc  cos  | 3jvg  (<*  -  *c  )  - 
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3c  ,  3,VC 
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sin  0c  Zc  sin  &NC  (<*  -  *c  )• 


(5.15) 


For  a  particular  d/lg ,  the  required  d  may  be  found  from  eq.  (5.15)  once  the  propagation  constants  /3c 
and  pNC  for  a  uniform  guide  have  been  found  from  Eqs.  (5.4M5.7). 


S3.  Results  and  Discussion 

Filter  designs  and  response  predictions  were  achieved  using  the  procedures  outlined  in  Section  5.2 
with  the  aid  of  computer  programs. 

A  scaled  image-guide  filter  was  built  and  tested  in  the  X-band  frequency  range.  The  predicted 
and  measured  insertion  lass  responses  and  the  dimensions  are  shown  in  Fig.  5.4.  Rexolite  1422  was  used 
as  the  dielectric,  due  to  its  favorable  electrical  and  mechanical  qualities.  To  enable  testing  with  a  net¬ 
work  analyzer  system,  a  suitable  transition  between  metal  waveguide  and  DWG  is  necessary.  The 
transition  used  consists  of  a  horn  and  an  H- plane  tapered  section  of  DWG.  This  has  been  shown  to  act 
as  a  low-loss  transition. 


The  predicted  response  is  for  matched  loads  at  midband.  (Maximum  ripple  across  the  band  due  to 
mismatch  is  less  than  1  dB).  Equation  (5.14)  predicts  the  lower  edge  of  the  stop-band,  so  a  slight 
adjustment  is  necessary  to  obtain  a  required  center  frequency.  The  measurements  agree  quite  well 
with  predictions.  Measured  data  in  the  11-12  GHz  range  deviated  from  the  predicted  response.  This 
deviation  was  due  to  radiation  in  this  frequency  range.  It  is  generally  accepted  that  grating-type 
slow-wave  DWG  structures  have  some  inherent  radiation  problems.  However,  the  radiation  region 
may  not  cause  a  problem  and  could  possibly  be  eliminated  by  tapering  the  steps  at  the  ends  of  the 
structure  [40].  Increasing  the  number  of  sections  or  the  step  size  increases  the  stop-band  insertion  loss. 
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Figure  5.4. 
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NORMALIZED  FREQUENCY  f/f0 


Predicted  and  measured  insertion  loss  for  grating  filter.  Rexolite  1422, 
hi  *  7  mm,ljfC/d  *  0.5,  d  =»  11.2  mm,  f  0  =  lOGffr,  w  *  20  mm. 

Predicted: 

- h  2  *  20  mm  \  N«20  sections 

- h  2  *  14  mm  J 

- h  2  =  14  mm  N“40  sections 

Measured: 

O  h  2  *  20  mm  > 

O  fc}*14  mm  y 


N-20  sections 


There  are  a  number  of  interesting  filter  structures  which  employ  gratings.  These  include: 

(i)  stagger-tuned  grating  sections  in  series  to  realize  a  low-pass  or  large  stop-band  band-reject  filter 
(Fig.  5.5(a)) 

(ii)  grating  with  coupler  (Fig.  5.5(b)).  The  forward  coupler  action  couples  the  reflected  power  from 
the  grating  to  achieve  a  band-pass  response  between  parts  1  and  2 

(iii)  tapering  in  the  grating  to  obtain  an  equal-ripple  response.  The  realizable  /fiNC  limits  this. 

A  simple  yet  effective  means  of  analyzing  dielectric  waveguide  grating  filters  has  been  presented. 
Experimentation  has  verified  the  theoretical  analysis.  This  technique  thus  becomes  a  useful  design  tool 
which  is  applicable  to  many  forms  of  periodic  and  non-periodic  guided  wave  structures. 
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Chapter  6:  CONCLUSION 

Several  millimeter-wave  waveguide  problems  have  been  investigated.  These  problems  involved  a 
single  or  a  series  of  discontinuities  in  fin-line  and  dielectric  waveguide. 

The  uniform  fin-line  was  analyzed  in  chapter  2.  The  problem  was  formulated  in  the  spectral 
domain,  and  a  numerical  solution  was  obtained  using  the  moment  method.  Results  for  the  propagation 
constant,  characteristic  impedance  and  mode  functions  were  given. 

A  two-dimensional  spectral  relationship  in  the  plane  of  the  metal  fin  was  used  to  formulate  the 
discontinuity  problem  in  chapter  3.  A  single  and  an  infinite  periodic  array  of  discontinuities  were  stu¬ 
died.  The  numerical  solution  was  obtained  by  using  a  moment  method.  Results  for  the  single  discon¬ 
tinuity  analysis  and  the  dispersion  curve  for  an  infinite  periodic  array  of  discontinuities  were  given.  In 
the  case  of  the  single  discontinuity,  efficient  numerical  evaluation  of  a  number  of  integrals  was  impor¬ 
tant.  Stability  of  the  solution  with  respect  to  variations  in  the  approximating  functions  was  a  prob¬ 
lem. 

In  chapter  3  the  fin-line  discontinuity  problem  was  formulated  with  unknowns  in  the  transverse 
plane  of  the  junction.  Two  numerical  procedures  were  used  to  solve  for  these  unknowns.  The  first,  the 
conjugate  gradient  method,  could  solve  the  uniform  problem  well,  but  would  not  converge  to  the 
correct  solution  in  the  case  of  the  discontinuity.  The  correct  solution  was  not  obtained  because  the 
number  of  degrees  of  freedom  in  the  unknown  was  greater  than  the  number  of  terms  possible  in  the 
Green’s  function.  The  second  approach  involved  a  use  of  the  generalized  variational  procedure,  which 
represented  the  unknown  in  terms  of  characteristic  functions.  Satisfactory  results  were  obtained  with 
this  approach. 

A  rectangular  dielectric  image  guide  grating  filter  was  analyzed  in  chapter  S.  The  effective  dielec¬ 
tric  constant  method  in  conjunction  with  a  transmission  line  model  was  used.  Satisfactory  agreement 
between  theory  and  experiment,  for  the  insertion  loss,  was  obtained. 
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